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I. INTRODUCTION 

Loop Quantum Gravity (LQG) [lj], a non-perturbative and background independent 
canonical quantisation of general relativity in four space-time dimensions, is one of the 
main approaches to quantising gravity. Even though the full theory of LQG is not yet com- 
plete, its successes lead us to apply LQG techniques in a simple setting where symmetry 
principles can be used. The application of LQG to the cosmological sector, known as Loop 
Quantum Cosmology (LQC) 0, has recently made significant progress. The difference 
between LQC and other approaches of quantum cosmology, is that the input is motivated 
by a full quantum gravity theory. The simplicity of the setting (typically homogeneity and 
then sometimes also isotropy, although recent progress has been made towards inhomoge- 
neous cosmologies jij), combined with the discreteness of spatial geometry provided by LQG, 
render feasible the overall study of LQC dynamics. 

Loop quantum cosmology is formulated in terms of SU(2) holonomies of the connec- 
tion and triads. In LQC, the quantum evolution is described by a second order difference 
equation, instead of the second order differential equation of the Wheeler-DeWitt (WDW) 
approach to quantum cosmology. As the scale factor increases, the universe eventually en- 
ters the semi-classical regime, and the WDW differential equation describes, to a very good 
approximation, the subsequent evolution. 

In the 'old' quantisation, the quantised holonomies were taken to be shift operators with a 
fixed magnitude, but later it was found @] that this leads to problematic instabilities in the 
continuum semi-classical limit. Indeed, as the universe expands, the Hamiltonian constraint 
operator creates new vertices of a lattice state, leading in LQC to a refinement of the discrete 
lattice. The effect of the lattice refinement has been modelled and the elimination of the 
instabilities in the continuum era has been explicitly shown 0]. Lattice refinement leads to 
new dynamical difference equations which, in general, do not have a uniform step-size. Thus, 
their study gets quite involved, particularly in the anisotropic cases. Recently, numerical 
techniques have been developed [1, [9} to address this issue. 

The correct refinement model should be given by the full LQG theory. One would have, 
in principle, to use the full Hamiltonian constraint and find the way that its action balances 
the creation of new vertices as the volume increases. Instead, phenomenological arguments 
have been used, where the choice of the lattice refinement is constrained by the form of the 
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matter Hamiltonian [7|. Here, we use another argument to specif, within a particular class 
(power law), the choice of the lattice refinement. 

There are many equivalent ways of writing the Hamiltonian constraint in terms of the 
triad and the holonomies of the connection, since at the classical level holonomies commute. 
However, each of these factor ordering choices leads to a different factor ordering of the 
WDW equation in the continuum limit. In what follows we explicitly demonstrate that the 
ambiguities at the classical limit of LQC, which is precisely the WDW equation, disappear 
only for a particular choice of a lattice refinement model. 



II. ELEMENTS OF LQC 

To quantise the gravitational Hamiltonian for isotropic flat cosmologies, we restrict our- 
selves to an elementary (fiducial) cell V, with finite fiducial volume; only in this volume 
spatial integrations will be performed. This is the usual approach to regularise divergences 
appearing in a quantisation scheme based on a Hamiltonian framework within flat homoge- 
neous models. 

Introducing a flat fiducial metric °q a b, in which the volume of V is Vo, the phase space 
variables p, c of loop quantum cosmology read 

H=vf^; c=X|_ T a, (2.1) 

(with the lapse function set to 1) where a is the cosmological scale factor and 7 is the 
Barbero-Immirzi parameter, labelling in-equivalent quantum theories. The triad compo- 
nent p, determining the physical volume of the fiducial cell, is connected to the connection 
component c, determining the physical edge length of the fiducial cell, through the Poisson 
bracket 

{c,P} = y, (2-2) 

(with k = 8nG) which (for this choice of variables) is independent of the volume factor Vq. 

The Hamiltonian formulation in the full LQG theory is based upon the Ashtekar vari- 
ables, namely the connection A\ and (density weighted) triad E®, arising from a canonical 
transformation of the ADM variables. Note that i refers to the Lie algebra index and a is a 
spatial index, with a, % — 1, 2, 3. They are given by 

-V 3 , ,i . T?a _ 



A\ = cV "^ w* ; E? = pVo ' V°? *i , (2-3) 

where °q is the determinant of the fiducial background metric °q a b = oj^ujui with uj l a a basis 
of left-invariant one-forms, and are the Bianchi I basis vectors Xf = Sf . 

After quantisation, states in the kinematical Hilbert space can be expressed as (linear 
combinations of) eigenstates of p, namely 

p|A*) = ^MI/*> , (2-4) 

which are diagonal, i.e. (/ii|/i 2 ) = 

Just as in full LQG, there is no operator corresponding to the connection, however the 
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action of its holonomy is well defined, 



ki\fi) = (Cs± + 2TiSn) \a) , (2.5) 



where Sn and Cs are given by 

1 / —ijlc ijlc \ — % 



Sn '^ = 2i ( 6 ^ ~ = T ^ + ^~ ^ ~~ ^ ' 

C^) = - (e^ + e^) - — + ft) + |/i - /i)) . (2.6) 



When /i is a constant (typically called _u is the 'old' quantisation), it is clear that the 
operator exp (iftc/2\ acts as a simple shift operator, namely 

exp (iftc/2^ = exp f^^J 1^) = I/" + • ( 2 -?) 

In the case of lattice refinement, ft — ft (//) is not a constant and the shift interpretation is 
no longer valid. However, one can change variables from fi to v: 

v = k lw)' (2 ' 8) 

for which _____ 

with k a constant. In these new variables the holonomies do act as simple shift operators, 
with parameter length k, for states labelled by u, defined as eigenvalues of / (p) (with / the 
implicit function giving /i (u); it is obtained by solving Eq. (12.81) ). Thus, 

e^\v) = e k ^\u) = \u + k) . (2.10) 

One should keep in mind, that the relationship between v and geometric quantities, such as 
volume, is more complicated than their relationship with ii. In what follows, we consider ft 
to be of the form ft = fi Q fi~ A , where _u is some constant In this case, one can explicitly 
solve Eq. (12. 8p to obtain 

/i (1 - A) 

This one-parameter family of lattice refinement models includes the 'old' (A = 0) and 'new' 
(A = —1/2) quantisations. Motivated by the full LQG theory, A is expected to lie in the 
range between —1/2 < A < [ 6j . There are several phenomenological and consistency 



arguments supporting A = —1/2 [lOj, [13J, nevertheless here we keep A as an undetermined 
parameter and show that the choice A = —1/2 leads to important consequences for the 
factor ordering of the continuum limit of the theory. 
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III. FACTOR ORDERING 



The classical gravitational part of the Hamiltonian constraint (with the lapse function 
set to unity) is given by 

1 f F! ai F! b i F k 

C grav = -- / d 3 x e ijk - n : : g , (3.1) 



r ./v 



Vldet^j ' 



where F„ 6 = d a A\ — d\,A\ + e] k A J a A^ is the curvature of the connection. Writing this in terms 



of our quantisable variables (p and the holonomies of c), it reads [11 

2i 



C 



grav 



K 2 foy 3 k 3 



(3.2) 



ijk 



There are many possible choices of factor ordering that could have been made at this point, 
since classically the action of the holonomies commute. Each of these possible choices will 
lead, in principle, to a different factor ordering of the resulting continuum WDW equation. 

We consider only factor orderings of the form of cyclic permutations of the holonomy and 
volume operators within the trace. They have the advantage of being trivially equivalent 
with respect to the spin indices, while they remain within the irreducible representations 
of the holonomies, hence avoiding the spurious, ill-behaved, solutions, present for higher 
representations 121 ] . Whilst these holonomies do not represent a complete set of factor 



ordering choices, nevertheless they include the choices commonly made in the literature 11 



12J . Finally, using these factor orderings we can demonstrate how different lattice refinement 
models alter the factor ordering of the resulting WDW equation in the continuum limit. 
Using Eqs. (I2.5p . fl2.6p . the action of the different factor ordering choices give us explicitly: 



e ijk tr (/ii/ij/ij 1 h 1 1 h k 



e ijk tr ( /ij/ij 1 h ] 1 h k 



e ijk tr (/i; l h s x h\ 



h[hj 



en k tr (h- 1 h 



eyktr ( h k 
eijktr 



h[h^h i h- /i k 



-2 lSii 2 Cs 2 ( CsV'Sii - SnV'CsJ 



24Sn Cs CsVSn - SnVCs Cs 



-12Sii 2 fCs^Sn - SnfCs") C S 2 



A2Cs ( CsVSn - SnVCs) Sn 



24Cs CsVSn - SnUCs Sn Cs 



-24 (CsVSn - Snt>Cs) SnW 
-24 [ Cs^/Sn - SnVCs ) Sn^ 2 



(3.3) 
(3.4) 

(3.5) 
(3.6) 
(3.7) 
(3.8) 
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where we have made extensive use of the trace identities: 

tr (n) 
tr (nrj) 
tr (Tirjrk) 

tr (TjTiTjTk) 



— tr (a-) — 
— tr Wj) 


n 








-tr (aio^) 




-1 

x 6jki 


1 , 

16 V J J 




-1 



1 



(3.9) 



Using Eqs. (I2.6P and defining V\u) = V v \u), the action of each of these factor orderings on 
the basis \u) can be calculated. For clarity, we will derive the continuum limit only for the 
factor ordering choice given in Eq. (13. 3p ; the other ones follow along similar lines and for 
completeness have been included in an appendix (Appendix A). 
Thus, 

Sn^Cs 2 (CsUS^ - SnVCs) \v) = ^ (V u+k - V v _ k ) [\v + Ak) - 2\v) + \u - 4Jfe)) . (3.10) 

Extending the above, one can obtain the action of the chosen factor ordering on a general 
state in the Hilbert space given by |^) = J2 u ipv\v)- The explicit difference equations for 
the coefficients ip u read 



eijktr ( /ij/ij/ij 1 h- 1 h 



I*) 



-3z 



E 



V v -zk — Vu-Bk ) ?Pv-4k — 2 I V u+ k — V v -k ) 



+(v u+5k - v„+3kji> v+Ak 



(3.11) 



We can now take the continuum limit of these expressions by expanding ip v « ip {y) as a 
Taylor expansion in small k/u, i.e. in the limit that the discreteness scale (k) is much smaller 
than the scale of the universe (which is given by v). By noting that the volume is given by 



KW~H^)] 3/ » 



(3.12) 



where /i (y) is obtained by Eq. (12.111) and we are neglecting a constant factor (n'yh/Q) 3 ^ 2 , we 
find 

3/[2(l-A)] 



K 



v±nk 



{y ± nk) a 



(3.13) 



where a — /Iq (1 — A) jk. 

In general, the above needs also to be expanded in the k/u — > limit and due to the k~ s 
factor in Eq. (13. 2p it is necessary to go to third order in both this and the Taylor expansion. 
To expand Eq. (13. lip we need the difference between the volume eigenvalues evaluated on 
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different lattice points, given by 

3k 



+mk 



2(1-A) 



a 3/[2(l-A)] l/ (l+2A)/[2(l-A)] 



1 — \~ 2^4. k / o o\ 

(n — m) + - r - r _ i — — [n — m j 



4(1 -A) v 



(1 + 2A) (4A - 1) /c 2 
24 (1 - A) 2 v 2 



(„• - m ») + o 



(3.14) 



Performing a Taylor expansion of Eq. (13.111) we get the large scale continuum limit of the 
Hamiltonian constraint: 



lim eijwtr ( h\h\h-. h, 

k/u^O J V J 



I*) 



-36i 
I - A 



a 



3/[2(l- 



,(1+2A)/[2(1-A)] 



d 2 ?/> 1 + 2A 1 d0 (1 + 2A) (AA - 1) 1 
d^ 2 + I- A ud^ + (I- A) 2 i^ 2 



(3.15) 



Taking A = reproduces the large scale factor ordering associated with the 'old' quantisa- 
tion, as expected 0]. Notice that A = —1/2, which corresponds to the 'new' quantisation, 
leads to the following very simple form of the evolution equation in the continuum limit: 



lim C grav |^) 

k/u-*0 



72/i //c f K n fh\ 



K 2 fry 3 



6 



d 2 ?/> 
d^ 2 



(3.16) 



where we have used a = 3//o/(2fc) and reintroduced all the constants. 

For fi = k we arrive at the following final result for the continuum limit of the WDW 
equation: 



lim Cgravl^) 



72 / K~fh\ 



k 2 ^7 3 



6 



\ 3/2 

) £ 



d 2 y> 
d^ 2 



(3.17) 



Repeating this tedious, but straight forward calculation for the other factor ordering 
choices given in Eqs. (13.41) - (13.8j) results in the following differential equations in the contin- 
uum limit (see the Appendix A for details), 



lim eHktr h\h: h 

k/u^O 1 y 1 



-36i 



hi) m 



x 



2A)/[2(1-A)] 



I- A 

d 2 ip 1 + 2A 1 d^ (1 + 2A) (44 - 1) 1 



-a 



3/[2(l-A)] 



dz/ 2 



+ 



1 - A z/ dz/ 



4z/ J 



W) , (3-18) 
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lim ejiktr ( h-, h, 



1-4 



a 



3/[2(l-A)] 



2A)/[2(1-A)] 



d 2 ^ 1 + 24 1 d?/> (1 + 2A) (44 - 1) 1 
d^" + 1-4 2^ch7 + (l - A) 2 



\v) , (3.19) 



lim eiiktr ( h, 1 ^ 
J V J 



lim e ijk tr h k 

k/u-*0 



lim r i|k tr( [k\v] kkk^K) |*I>) 



-36z 



1-4 



-a: 



.! (2(1-A)) fc 3 ^(1+2A)/(2(1-A)) ^ | v 



(3.20) 



Once again we see that 4 = —1/2 results in a particular simplification, in that all of 
the considered factor ordering choices reduce to Eq. (I3.17p . This is to be expected since 
quantum factor ordering ambiguities should disappear at the classical level. The crucial 
finding is that the LQC ambiguities, associated with the factor ordering in the Hamiltonian 
constraint disappear at the continuum described by the WDW equation, only for a lattice 
refinement power law model with 4 = —1/2. Thus, only this model has a non-ambiguous 
continuum limit. 

In addition to this lattice refinement model providing a unique choice of factor ordering 
for the continuum limit equations, the action of the volume operator is greatly simplified, 



Vu+nk - V u+mk = (n - m) 



3/i c 



(3.21) 



This is no accident since the initial motivation for this quantisation procedure was that 
the volume, rather than the area, should get quantised. The consequence is that quantum 
corrections to this classical equation (i.e. quantum corrections to general relativity) enter 
only in the Taylor approximation i\){y + nk) m ip{v) + . . . and not in the expansion of the 
volume terms, i.e. Eq. f!3.14p requires no approximation. 

Finally, to relate this result to more usual variables, we can use /x ~ p = a 2 and v ~ /i 3 ^ 2 , 
to find that the factor ordering of the Wheeler-DeWitt equation predicted by the large scale 
limit of LQC reads 

A2J, a / aj.\ 

(3.22) 



Cgrav ~ 



dis 2 da V da 



where constants, but no factors on a, have been dropped. 
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IV. CONCLUSIONS 

Just as in the quantisation of standard fields, loop quantum cosmology results in factor 
ordering ambiguities. One expects that the classical limit should unambiguously result to 
the original classical equation. Here, we have shown that in general this is not true at the 
level of the Wheeler-DeWitt equation, which is the classical limit of LQC (it is obtained 
when the discreteness scale is set to zero). This ordering ambiguity disappears however for 
the particular lattice refinement model given by A = —1/2, which is typically called 'new' 
or 'improved' quantisation. 

The work presented here can be viewed in two ways: One could accept the phenomeno- 
logical 0, 10] and consistency 13] requirements indicating that A = —1/2 is the correct 



quantisation approach. In this way, we have shown that LQC predicts a unique factor or- 
dering of the Wheeler-DeWitt equation in its continuum limit (at least for the particular 
class of factor ordering considered here). Alternatively, one could require that factor or- 
dering ambiguities in LQC should disappear at the level of the Wheeler-DeWitt equation. 
In this case, we have shown that the lattice refinement model should be A = —1/2. In 
either case, we have clearly demonstrated that there is a strong link between factor ordering 
and lattice refinement in LQC; these two ambiguities are closely related. In fact, we have 
shown that by specifying a particular lattice refinement model we can uniquely determine 
the factor ordering of the equation in the continuum limit, and vice versa. 

In conclusion, it is remarkable that the requirement for the Wheeler-DeWitt factor or- 
dering to be unique, is precisely the same requirement reached by physical considerations of 
large scale physics and consistency of the quantisation structure. In particular, it has been 
previously shown [l(| that for LQC to generically support inflation and other matter fields 
without the onset of large scale quantum gravity corrections, A should be equal to —1/2. It 
has been recently shown [l3[ that physical quantities depend on the choice of the elementary 
cell used to regulate the spatial integrations, unless one chooses A = —1/2, and that sensible 
effective equations exist only for this choice. Taking this together with the uniqueness of 
the factor ordering of the Wheeler-DeWitt equation, which as we have shown also requires 
A = —1/2, it is clear that several vastly different unrelated approaches have converged on 
the same restriction of the theory. It is possible that this restriction can be used to improve 
our understanding of the underlying full loop quantum gravity theory. 
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APPENDIX A 



The action of the different factor ordering possibilities considered in Eqs. (I3.4l) - fl3.8p on 
a basis state \v) is given by 



S^Cs (CsFSn - SnFCs) Cs\v) = -J 



Sli 2 (CsFSn - SnVCs) CsV) 



(V v+2k ~ V v ) (\v + 4k) -\u + 2k) - \v) + \u - 2k) 



—% 
32 



+2 (y u+k - V u - k ) (\u + 2k) - 2\v) + \v - 2k) 
+ (K-fc - K- 3 fc) (V) - 2\u - 2k) + \u- 4k) 



Cs 2 (CsVSn - SnVCs) SnV) = ^ 



(V v+3k - V u+k ) (\v + 4k) + 2\v + 2k) + \v) 



-2 (K+fc - V v - k ) (\v + 2k) + 2\v) + \v - 2k)^j 
+ (K-fc - K- 3 fc) (V) + 2\u - 2k) + 4k)} 



+ (K - K-2fc) (\u + 2k) - 1 1/) -\u- 2k) + \u- 4k)) 

(Al) 

{V v+3k - V v+k ) (\u + 4k)-2\u + 2k) + \u) 



(A2) 



(A3) 
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C s ( CsFSn - SnFCs) Sn 2 Cs|i/) = -J 



(K+4fe - V v+2k ) (V + 4fc > + W + 2k )) 

-(K +2fe -K) {\u + 2k) + \u)^j 
-(V v -V v -2k){W) + W-2k)"j 
+ (K- 2fe - V„-4k) (\v-2h) + \v- 4A;)) 



(A4) 



Csrsn - Sure* ) Sn 2 CsV) = ^ 



(A5) 



Extending these to a general state in the Hilbert space given by |\&) = X^Vvl* 7 ) gives, 
e ijk tr (hfi^h^hk /i k \ V h^j = ^ ^ f K-2fc - K-4fc) ^„-4* - (K - V v - 2 kj^u-2k 

-^K+2fe - + (V v+4k - V u+2 k)lpv+2k 

+ (V v -2k ~ K-4fc)?/V-2fc - (V v ~ V v -2k)4> V 

-(y v+2 k - v^)%i) v +2k + (y v+ 4k - v v+ 2kjip I4A6) 



e ijk tr (h^hfh. h^,V kkj\^) 



-3i 



-4^K+fc - V v -k)il) v + (y v - k - K- 3 fc)^ 

+ (y u+3k - K+fc)^+4fcj , (A7) 
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e ijk tr (h7% [h k \ kkK 1 ) l*> = -j- Yl ( ( V » ~ V v-2k)*Pu-4k + (V u+2k - vjjip v . 



2k 



-3i 



V u - V v -2k)lpv-2k ~ (V u+2 k ~ V u ^lp u 

V v - K-2fcJ^ - (v u+2 k - V„)lf) 



u+2k 



+ (v u - v v -2k)i>v+2k + (v u+2 k - Vvjifjv+Ak ) kXAs) 



eijktr ( h k 
eyktr 



Performing the Taylor expansion of these as in Eq. (13.151) we get, 



(A9) 



lim e;jktr ( h,h- hi 

k/u-,0 \ 3 



-36z 



V 



(1+2A)/(2(1-A)) 



I - A 

d 2 ^ 1 + 2A 1 ch/> (1 + 2A) (44 - 1) 1 



a 



3/(2(l-A)) 



dz/ 2 1 — A v dv 



(1 - A) 2 4z/ 2 



W , (Aio) 



lim eijktr ( h: h: 

k/u^O J y J 



36i 



x 



* 3 £ 



(1+2A)/(2(1-A)) 



1 - A 

d 2 ^ 1 + 1 dip (1 + 2A) (44 - 1) 1 



-a 



3/(2(l-A)) 



dz/ 2 + I — A 2v dv 



+ 



(1 - A) 8z/ 



i>{v) \\v) ,(A11) 



k/ 

1: 

k/ 



k/ 
-36i 



lim e ijk tr (fy ^ k ft,^ 1 , V hjijhl ^ |^) = 



-a 



3/(2(1- 



A 



(1+2A)/(2(1-A)) 



d 2 ^ 



